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A system of equations governing the structure of a steady, relativistic radiation dominated shock 
is derived, starting from the general form of the transfer equation obeyed by the photon distribution 
function. Closure is obtained by truncating the system of moment equations at some order. The 
anisotropy of the photon distribution function inside the shock is shown to increase with increasing 
shock velocity, approaching nearly perfect beaming at upstream Lorentz factors r_ >> 1. Solutions 
of the shock equations are presented for some range of upstream conditions. These solutions are 
shown to converge as the truncation order is increased. 
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Relativistic shocks play an important role in essentially 
all classes of high-energy compact astrophysical systems. 
Behind these shocks the bulk energy of a coUimated out- 
flow expelled from the central engine is dissipated and 
converted to the radio-through-7-ray emission observed. 
The micro-physics of the shock depends on the conditions 
in the upstream and downstream regions. Shocks that 
form in a region of small optical depth for Thomson scat- 
tering are mediated by collective plasma processes, and 
are termed coUisionless shocks [l|. Various observations 
suggest that these shocks may provide sites for accelera- 
tion of nonthermal particles (both electrons and protons) 
and, at least in some systems (e.g., GRB afterglow), may 
somehow generate macroscopic, sub-equipartition mag- 
netic fields in the downstream region [12]. The latter 
processes may in fact be an inherent part of the shock 
structure Q. The micro-physics of coUisionless shocks 
is highly involved. Recent efforts have led to important 
progress in the study of coUisionless shocks, e.g., Refs. 
0, S H H 0| ' but despite these efforts many key issues 
remain as yet unresolved. 

When the ratio of radiation to matter pressure in the 
post shock region exceeds a certain level the shock tran- 
sition may become mediated by photons. Such condi- 
tions are anticipated in, e.g., internal 8, 9] and oblique 
jlQ] shocks producing the GRB prompt emission, during 
shock breakout in supernovae and hypernovae [111 Il2| . 
in blazars, and in accretion flows onto a black hole 
[13-, 14, 15|]. The effect of radiation on the structure of a 
non-relativistic, strong shock has been studied by differ- 
ent authors [ll, [13, [M E H HH, HI and appUed to var- 
ious astrophysical systems, such as supernovae \il\ and 
accretion flows onto neutron stars [i8|. The structure 
and transmitted photon spectrum of a cold, radiation 
mediated shock was computed in Ref. [2^, where it was 
shown that bulk Comptonization on the converging flow 
can give rise to a power law spectrum at high energies 
with a spectral index that tends to unity for large Mach 
numbers. The model has been generalized later to incor- 
porate thermal effects [21, 22|. The analysis in Ref [20] is 
restricted to sufficiently optically thick shocks with small 
upstream velocity, viz., /3_ << 1. In this regime the 



scattered photon distribution is nearly isotropic at every 
point inside the shock and the diffusion approximation 
can be employed to solve the transport equation describ- 
ing the evolution of the photon population across the 
shock [13] ■ Moreover, to order /3^ the radiation field sat- 
isfies the equation of state Prad = Urad/'^, which provides 
a closure condition for the set of hydrodynamic equations 
governing the shock structure. This simplifies the anal- 
ysis considerably, but renders its applicability to most 
high-energy compact sources of little relevance. 

Several complications arise in the relativistic regime. 
Firstly, the optical depth across the shock, r ~ c//3_, 
approaches unity as the shock becomes relativistic, hence 
photons do not experience multiple scattering. This is a 
consequence of the fact that the average energy a photon 
gain in a single scattering is large, viz., Ae/e > 1. As a 
result, the photon distribution function across the shock 
is anticipated to be highly anisotropic, rendering the dif- 
fusion approximation inapplicable. Obtaining a closure 
of the hydrodynamic shock equations then becomes an 
involved problem. Secondly, the optical depth of a fluid 
slab having a Lorentz factor F > 1 depends on the angle 
9 between the photon direction and the shock velocity as 
dr = F(l— /3 cos 9)dx, and needs to be properly accounted 
for. Thirdly, pair creation may become important if the 
photon energy exceeds the pair creation threshold. Suf- 
ficiently relativistic shocks may therefore become domi- 
nated by electron-positron pairs. From the physical point 
of view, then, a relativistic radiation-dominated shock is 
a complex physical phenomenon, posing an intriguing un- 
solved problem. In what follows we derive a coupled set 
of hydrodynamic equations that provide a complete de- 
scription of the shock structure, starting from the general 
form of the transfer equation obeyed by the photon dis- 
tribution function. We then propose a method to close 
the resultant system of equations and present solutions 
for some range of upstream conditions. Convergence of 
these solutions is verified. 

The fluid in the shock transition layer is a mixture 
of baryons, e^ pairs, and radiation. We denote by 
u°' = (F,F/3) the 4-velocity of the mixed fluid, by nt 
and n± the number density of baryons and pairs, respec- 



tively (the total number of electrons in this notation is 
Tie = nb + n_), and by T^^", T^", and T/"" the stress- 
energy tensors of baryons pairs and radiation, respec- 
tively. The entire system must conserve energy, momen- 
tum and baryon number: 
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The energy-momentum tensor of the radiation can be 
expressed explicitly in terms of the photon distribution 
function, fr{k,x'^), as 



T/'-(a;^)= / k^'k\fr{k,x^') 
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where /c^ denotes the 4-momentum of a photon. The 
distribution function satisfies a transfer equation, which 
we express in the form [23| 
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Here n'^ — n'j^ + n'_^_ + n'_ is the net number density of 
scatterers (electrons plus positrons), as measured in the 
fluid rest frame (henceforth, primed quantities refer to 
the fluid rest frame), R{k,ki) is the redistribution func- 
tion for Compton scattering from initial state ki to a 
final state k, the operator Cpp{fr, f±, k) accounts for the 
change in fr due to e pair creation and annihilation, and 
Sk is a source term associated with all other processes 
that create or destroy photons (e.g., free- free emission 
and absorption) . Conservation of energy, momentum and 
number of quanta of the interacting pair-photon system 
implies 

d^jr = -Jk'^Cpp{frJ±,k) — , (5) 

g^Ku'^)-~JCpp{frJ±,k)^. (6) 



The above set of equations augmented by appropriate 
boundary conditions upstream provides a complete de- 
scription of the shock transition layer. To simplify the 
problem we ignore KN effects in eq. (|4]). We then have 

Rik',k'^)^^[i+{k'.^r]s{k'"-kr), (7) 

where ctt is the cross section for Thomson scattering. 
Multiplying Eq. Q by fc'^, integrating over final states k, 
using Eq. ([7]) , and relating k' and k through the Lorentz 
transformation k'^ — hKk'^, we obtain 
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We consider a plane-parallel shock in the (y, z) plane 
(the flow is moving in x-direction), and solve the above 
equations in the shock frame, where the system is as- 
sumed to be in a steady state. For the range of up- 
stream conditions explored here the downstream tem- 
perature is well below WeC^, and so the pairs are sub- 
relativistic. To a good approximation we then have 
T!^ = n±meC^uf^u'^ . Since the shock is radiation dom- 
inated we can neglect the pressure contributed by the 
particles to get T^^" = mpC^n'i^u'^u'^ . Using the latter re- 
sult and the transformation T/o^ = A°(-/3)A^(-/3)T'^'', 
the zeroth component of Eq. ^ yields an equation of 
motion for the fluid Lorentz factor: 
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where dr = cryfijEdx is the angle-averaged optical depth 
for Thomson scattering of a fluid slab of thickness dx, 
and J = m.pC^n'f^u'' is the integral of Eq. ((T|) . To com- 
plete the set of equations a closure relation is needed. 
Far upstream {x — — oo) and downstream {x = +oo) the 
radiation field is in equilibrium, and the equation of state 
•yj^ixx _ rpiQQ applies. However, within the shock transi- 
tion layer this relation no longer holds. It is convenient 
to define a dimensionless function ^ that measures the 
deviation from isotropy at any given point: 



T;-- = i(i + Or; 
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In general ^ is a function of the shock velocity and its 




FIG. 1: ^c versus upstream velocity /3_, for different values of 
the ratio of radiation pressure and rest mass energy density, 
as indicated 

derivatives, that is, ^ — ^(u^,-u^). Owing to relativis- 
tic boosting the radiation field is expected to be beamed 
preferentially in the direction of fluid motion. Conse- 
quently we anticipate values between ^ = (complete 
isotropy) and ^ = 2 (perfect beaming in the x direction) 



at any given point in the shock transition layer. Let us 
denote by Ci and C2 the integrals of the zeroth and x 
components of Eq. ([2]), respectively. These integration 
constants are determined by the conditions upstream the 
shock, at a; = —00. From Eqs. ([1]), ([2]) and (fTO|) we ob- 
tain the net photon flux in the local rest frame of the 
fluid: 
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It can be readily shown that for ^ = the flux T/"^ has a 
root at /3 = /?+ , and f] — (3-, where /?+ (/3_ ) is the velocity 
of the downstream fluid, as determined from the shock 
jump conditions. It is also seen that the comoving photon 
flux has a singular point at /^^ = (1 + ^)/3. Since T^P^ 
must be finite everywhere inside the shock it implies that 
the numerator on the R.H.S of Eq. pT|) must vanish at 
the critical point. For upstream velocity /3_ > l/\/3 we 
can then readily constrain S,- For illustration, we present 
in Fig. 1 plots of £,c, the value of ^ at the singularity, for 
cases where the energy density of pairs can be neglected 
in Eq. (fTTj) . The presence of pairs is not expected to 
change the results significantly. The maximum value of 
^ inside the shock must lie between ^c and 2. As seen 
from Fig. 1, ^ approaches the maximal value ^ = 2 as 
the shock becomes relativistic, implying a highly beamed 
radiation field inside the shock. 

In the limit /?_ < < 1 we expand F in powers of f3 and 
T/°^ in powers of /3 and £,. By employing Eqs. pT|) and 
^ and noting that the average photon energy is well 
below the pair production threshold, in which case we 
can set n± = 0, we recover, to order 0{P'^), the result 
derived in Ref. [M: 
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where /i = PI P- is the normalized fluid velocity and 
7r_ — {T^^/J(3)^ is the ratio of radiation pressure and 
ram pressure far upstream. It can also be readily shown 
that to this order T°^ = dT^^/dr, implying that in this 
limit photon transport is indeed a diffusion process. 

In the relativistic case the particle content of the 
shocked fluid may be dominated by e pairs. The pair 
density can be computed from Eq. ^ once the collision 
term Cpp is known. However, the latter depends on the 
energy distribution of photon, the calculation of which is 
beyond the scope of this paper. As a rough estimate of 
the pair density we use the equilibrium value. 
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To complete our treatment we need to derive an equation 
describing the change in ^ across the shock. We do so 
by integrating Eq. ^ over fc" and (j>. Defining /(/i) = 



/ {k'^ffdk°d4>, and likewise /'(/x') = / {k'°)^f'dk'°d(t)', 
recalling that I{ii) — r'*(l + f]fi')'^I'{^'), and noting that 
dfi'/dr = {fi' - l)(r/3)"Mr/dT, we arrive at. 
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where n'larKpp = r3(l - f]^)^ J {k")^Cppdk"d(j). We 
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FIG. 2: Fluid velocity, comoving photon flux T',?^ , and 
anisotropy parameter ^ as functions of optical depth, for mat- 
ter dominated (left panel) and radiation dominated (right 
panel) fluid far upstream. The blue lines in both panels are 
solutions obtained to second order and the red lines to third 
order (see text for further details) 

solve Eq. (|14p by expanding /' and Kpp in terms of 
Lagendre polynomials: 
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/'(a^') = -I]r,fePfc(M'); ^vvkiJ^) = 2^Cfe^fc(M')- (15) 

Note that with this normalization T/"" = 7?o, T;°^ = 771/3 
and ^ = 2772/5770- Substituting Eq. ^^ into Eq. (|Ti)) 
yields an infinite set of ODE's for the unknown variables 
r\n- To simplify the analysis we keep only the terms Co 
and Ci in the expansion for Kpp which, using Eq. ([5]), 
can be expressed as Co = T'^^PdT^jdT - dT'^/dr) and 
Ci =F2(/3dT^^/dT-dr|^/dr)[2i]. A closure relation can 
be obtained by truncating the expansion of /'(/^') in Eq. 
(fT5|) at some order n. It can be shown that the resul- 
tant set of equations has singularities at certain values 
of /?, that are independent of the conditions upstream. 
The root in the denominator of Eq. pT|) , which can al- 
ternatively be derived using the moment equations, is 
one of these points. As seen, the full transfer equation, 
Eq. (fT4|) . has a singularity only at the point /3 = —11' 
(corresponding to /z = in the shock frame), where the 
optical depth tends to infinity. This implies that any 
physical shock solution of Eq. P^ must automatically 




FIG. 3; Velocity profile for an upstream fluid with r_ — 2 
and different values of Tr'^'^/J. 



pass through the singular points of the moment equa- 
tions. This is demonstrated already in Eq. (jlip where 
the solution for ^ in all cases studied automatically takes 
the value at the singularity that leaves T/^^ finite. These 
singular points are therefore different in some respects 
from the critical points of a hydrodynamic system that fix 
physical parameters of a transonic solution (see Ref. [241 
for a detailed account). We were able to find physi- 
cal shock solutions of the moment equations that pass 
smoothly through these singular points, for mildly rela- 
tivistic upstream velocities (r_ ^ a few). Full details will 
be given elsewhere. For the cases studied below we find 
that the solutions converge as the order of truncation is 
increased. An example is shown in Fig. 2, where the fluid 
velocity /3, comoving photon flux T^°^, and anisotropy pa- 
rameter ^ are plotted as functions of the angle averaged 
optical depth, r, for upstream velocity f3- = 0.7. Second 
order (blue lines) and third order (red lines) solutions 
are compared in each panel. As seen, for the radiation 
dominated upstream condition (right panel) the agree- 
ment is perfect. For the matter dominated condition the 
convergence is slower. This is in part a consequence of 
the larger anisotropy (larger values of £, and T^°^) in the 
latter example. We have obtained also higher order solu- 
tions for this case and found nearly complete convergence 
at fourth order. Shock profiles computed for upstream 
Lorentz factor r_ = 2 are exhibited in fig 3. The ultra- 
relativistic regime requires analysis of higher order terms 
and is currently under study. 

From the above analysis it is evident that an appre- 
ciable fraction of the shock energy should be converted, 
via bulk Comptonization, to high energy photons having 
energies well in excess of the thermal peak. If the optical 



depth in the upstream region is not too large then these 
photons would escape before being thermalized. Con- 
sequently, a non-thermal spectral component appears to 
be an inherent feature of relativistic radiation dominated 
shocks and requires no particle acceleration, as in the case 
of collisionless shocks. The transmitted spectrum should 
extend up to an energy of ~ rUeC^, as measured in the 
shock frame, above which it will be suppressed by KN 
effects. In the frame of the observer the high energy cut- 
off of the emitted spectrum will be boosted by a factor 
Ts, the shock Lorentz factor. For GRBs this implies cut- 
off energy that can exceed 50 MeV or so in the observer 
frame. This mechanism can therefore account for the 
spectra observed in most GRBs. Detailed calculations 
of the transmitted spectrum are underway (Bromberg & 
Levinson, in preparation). 

Relativistic photon mediated shocks can provide a 
means for producing nonthermal spectra also in purely 
leptonic fireballs, the consideration of which is motivated 
by recent post-SWIFT discoveries of a shallow afterglow 
phase at early times [25[. According to some interpre- 
tations (e.g., Ref. [2y]) these observations indicate pro- 
longed activity of the central engine that, if true, implies 
that 7-rays are emitted during the prompt phase with 
very high efficiency. Such episodes can be most naturally 
explained as resulting from photon mediated shocks in a 
pure electron-positron plasma. Further discussions can 
be found in Refs. jlO . \2T \ 

As seen from Figs 2 the scale of a photon mediated 
shock is a few Thomson depths, which for parameters 
typical to most compact, relativistic systems is several 
orders of magnitudes larger than any microphysical scale 
(e.g., the skin depth and Larmor radius) associated with 
collisionless shocks. Thus, unlike collisionless shocks, 
photon mediated shocks cannot accelerate particles to 
nonthermal energies. This has direct implications for 
production of cosmic rays and VHE neutrinos in those 
systems. A particular example is the recent proposal that 
failed GRBs may be prodigious sources of TeV neutrinos 
[28| . The idea is based on the ad hoc assumption that 
internal collisionless shocks that form in the chocked out- 
flow accelerate protons to very high energies. However, 
the large Thomson depth anticipated in the region where 
the shocks form |28i] should render them radiation dom- 
inated, as our analysis indicates, and so they may not 
be able to provide the required sites for the acceleration 
of protons. Thus, the effects of radiation domination are 
likely to have a strong influence on the predictions pre- 
sented in dll. 
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